APPLICATIONS REFERENCE LIBRARY; 

HYBRID COMPUTATION 



HYBRID COMPUTER ANALYSIS OF ELECTROCARDIOGRAPHIC DATA 



INTRODUCTION 

This Study describes the application of the EAI 
HYDAC® 2000 (HYbrid Digital-Analog Computer) 
System to the analysis of electrocardiagram (EKG) 
data. The type of analysis employed in the study 
represents one of a class of data reduction tech- 
niques which can be instrumented effectively on this 
general pur^DOse hybrid computer. 

The overall area of data analysis is an essential 
one to the Bio-Medical Engineeringprofession since 
most of the work being done at present is on an 
empirical basis. Thus, large quantities of data, 
representing measurements from various sections 
of the lx)dy, must be gathered, and detailed analyses 
of this data, as is necessary to establish significant 
correlation between physiological malfunctions and 
data irregularities, must be performed. 

A typical but self-limiting approach to a project of 
this type has been to make use of special purpose 
equipment designed specifically for cardiogram 
analysis, for autoradiographic analysis, for pul- 
monary circulation transport studies or, indeed, for 
investigating any of the physiological processes 
within the scope of bio-engineering. The basic 
limitations of this approach, of course, are: 1) 
the impossibility, or very great difficulty, of 
adapting the equipment for any data analysis other 
than that for which it is designed, and 2) the 
extensive hardware requirement for obtaining val- 
uable but, at best, partial information. 

The HYDAC program to be described (block dia- 
gram shown in Figure 1) is an economically-attrac- 
tive alternative approach to satisfying the need for 
physiological data reduction in general, and electro 
cardiographic data analysis in particular. Further- 
more, since the HYDAC computer is a general pur- 
pose computer, there is the additional capability 
for complete system simulation of physiological 
functions. As in industry, where the mathematical 
^ ' simulation of processes and mechanisms is rapidly 
being adopted as a standard procedure in moderni- 
zation and optimization projects, so in bio-en- 
gineering is such system model-building assuming 
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a place of importance. In this way, the computer 
can provide still greater insight into the operation 
of the human body by making possible, say, the 
analysis of subsystem interaction in the function 
of the entire organism or the introduction of some 
control functions. 

Indeed, these two aspects of computer operation — 
data analysis and direct simulation — enhance each 
other in the furtherance of bio-engineering accom- 
plishment. That is, the development of mathematical 
models for simulation will require the analysis of 
vast quantities of empirical data. Analysis of the 
model, however, will point out areas about which 
more or different information is needed in order to 
complete the mathematical description. Thus, the 
quality of data would be upgraded and the mathe- 
matical model continuously refined to the point 
where the simulated system most closely ap- 
proached the actual physiological process. 

PROBLEM STATEMENT 

There is a great deal of current interest in methods 
for analyzing human electrocardiograms. The basic 
signal to be analyzed (e.g., a magnetic tape record) 
is a measure of the voltages determined by three 
or more leads placed in contact with the skin of a 
human body. These voltages are generated within 
the heart, and the specific sequence of the pulses 
generated is indicative of stages in the expansion- 
contraction cycle of heart action. Actually, the 
sequence is produced by an electrical polarization 
and subsequent depolarization of the muscle tis- 
sue which contracts during the pumping action of the 
heart. The height of the pulses so generated and the 
time intervals between them are affected by various 
disease states of the human body. In order to de- 
rive statistical correlation between these pulse 
measurements and human disease — in order, that 
is, to make the cardiographic analysis meaningful- - 
it is necessary that a large number of cardiograms 
be analyzed in minimal time so that the information 
obtained corresponds as closely as possible with 
actual, real-time heart action. This is of special 
importance in physiological studies where the 
cardiographic analyses are part of an on-line 
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PROGRAM n: SIMPLIFIED DIAGRAM.PARAMETER MEASUREMENT 
Figure 1: Block Diagrams of Programs I and II of HYDAC Cardiographic Data Processor 



system in a complete system simulation, and must 
proceed in real-time to keeppace with complemen- 
tary and/or supporting bodily functions. Further- 
more, in order to remove possible bias in the 
analysis of cardiograms as might be desireable, 
say, in insurance investigations where the bias 
of advanced age in the subject, for example, might 
cause an over-critical interpretation of a trace, it 
is necessary to standardize the measurement of 
the critical pulse heights (amplitudes) and time dura- 
tions in the recorded waveform. 



The HYDAC program described here was pre- 
pared to satisfy these two requirements, viz: to 
process cardiogram records at essentially the 
same rate that they are generated by human pa- 
tients, and to establish criteria by which critical 
pulse heights and time durations can be measured 
and interpreted consistently. The program is divid- 
ed into two sections: Program I preprocesses the 
data in order to circumvent a basic "noise" 
problem — signal distortion caused by spurious EKG 
baseline variations; Program II performs the 
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actual pulse measurements. Figure 1, above, shows 
a simplified block diagram of the HYDAC Cardio- 
graphic Data Processor. 

The diagram makes evident one additional benefit 
of this general purpose hybrid approach. The use of 
two separate programs has the advantage that the 
data smoothing operation in Program I produces 
its own paper tape output which can be processed 
further in a variety of ways not considered by 
Program II, thus adding to the versatility of the 
entire computer operation. 

PROGRAM I: Data Pre-Processing 

A. Signal Composition — generation of a standard wave- 
form: The various systems for locating leads on 
the human body for EKG examinations make possible 
a number of different methods for obtaining a com- 
posite signal from the several voltages measured. 
The objective, in any event, is to relate this com- 
posite signal to the ''resultant heart potential 
vector"... the quantity involving direction as well as 
magnitude of the average polarization potential of 
the heart as a function of time. 

For illustrative purposes, an orthogonal lead sys- 
tem will be considered in this paper. For an 
orthogonal lead system, the leads are placed in 
positions on the body so that the resulting amplitude 
of the heart vector is obtained by the simple 
formula 

EKG3D =>/x2 + y2 + z2 (1) 

where EKG3D = amplitude of resultant waveform 
X, y and z = amplitudes of voltages obtained 
from orthogonal lead system 

This formula destroys information about the phase 
(or sign) of the voltages x, y and z and, since the 
phase inversion of a particular pulse is a significant 
medical event, some method of remembering the 
phase must be introduced. For example, either a 
logic signal indicating plus (+) or minus (-) for 
each X, y and z wave must be produced, or else 
(as done in this study) x, y and z can be shifted by 
a small amount until they become voltages which 
are always positive. 



By using the shifted waveforms, i.e., x + ej^, 
y + e2, and z + e3, the phase of the resultant will 
follow the dominant phase of x, y and z at the ex- 
pense, only, of a small amount of signal distortion. 
The process is performed by an analog computer 
producing 



EKG3D =V (X + e^)^ + (y + e^)^ + (z + e^)^ (2) 

which, for specific standard valuesof ei, e2 and e3, 
is the basic single waveform processed by the com- 
puter. 

B. Data Smoothing: A primary problem in the ac- 
curate determination of pulse heights and dura- 
tions is the elimination of spurious baseline 
variations in the electrocardiogram. 
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Figure 2: Typical EKG Baseline Variations 

Since it is the height of the pulse above the variable 
baseline which is desired (Figure 2), it is necessary 
to remove the effect of baseline eccentricities. 
Simple filter circuits fail to remove this baseline 
component without altering significantly the wave- 
shapes to be measured. Consequently, an averaging 
scheme was adopted instead, in an attempt to re- 
move these variations without distortion. This task 
is the primary one performed by Program I, 

The basic scheme, as shown in Figure 3, is to 
operate on the almost-periodic waveshape of the 
EKG* after resetting the zero reference level 
after each period. This resetting procedure re- 
moves a considerable portion of the long term 
baseline variations. However, short term varia- 
tions exist within periods and, to remove these, 
the HYDAC program averages each period of 
the waveform with all previous wave periods. 

In this manner, the remaining baseline variations 
are averaged out to produce a clear, noiseless, 

* Actually, 3 cardiogram potentials have been combined to form 
EKG3D as described previously. 
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Figure 3: Basic Averaging Scheme Used in Program I 

average representation of a cardiogram waveform 
which is punched out on paper tape. A large 
number of recoixis can beprocessed in this manner. 
The resulting averaged waveforms then can be 
inputted either directly to Program 11 to complete 
the EKG analysis or to other devices for further 
processing, as desired. 

PROGRAM II: Parameter Measurement 

Program II accepts the paper tape output of Pro- 
gram I and stores it in 'memory'. The average 
waveform then is played back through a digital- 
to-analog P/A) converter into an analog com- 
puter where the occurence of the various pulses 
is detected, and their amplitudes and appropriate 
time intervals are measured. An analogdifferentia- 
tor detects the peaks and creates signals controlling 
time-interval measurement by integrators. 

This data is transferred to the Digital Operations 
System (DOS) of the HYDAC Computer and again 
punched onto paper tape. A series of logic controls 
in the DOS oversees the measurement process on 
the analog computer, detects the need to abandon 
the search for a particular pulse (i.e., determines 
that this pulse is not present), and begins expecting 
a succeeding type of peak. These non-detected peaks 
are remembered and an appropriate tape code in- 
dicates that their measurement has not been per- 
formed. 




ABANDON POINTS 
A I : ABANDON SEARCH FOR P PULSE 
A 2*. ABANDON SEARCH FOR Q PULSE 
A3: ABANDON SEARCH FOR S PULSE 
A4: ABANDON SEARCH FOR T PULSE 

Figure 4: Pulse Search Procedure 



Measurement is made of five pulse amplitudes, 
five basic time periods and, where detected, the 
variance of critical time measurements (defined 
below). These analog-calculated quantities then 
are punched onto a paper tape output record along 
with appropriate patient identification. An analog 
plot of the average waveform is generated also as 
a further piece of output data. 

GENERAL PROGRAMMING CHARACTERISTICS 

Certain characteristics of the EKG record, as well 
as a desire for real-time processing, dictate spe- 
cial programming restrictions in the parameter 
measurement phase of the program. The most basic 
characteristics of the cardiogram are that it con- 
tains five separate pulses labelled P, Q, R, S and T, 
and that only one of these (i.e. the R pulse) is 
necessarily present and readily detectable inhuman 
cardiograms (Figure 5). The P, Q, S and T pulses 
may be extremely small or not present at all in 
some cases. Thus, schemes which detect the R 
pulse are the only ones considered reliable on the 
unaveraged waveform. This means that a given 
point on the R pulse (say, a point of given negative 
slope) must be used as a reference and that all 
averaging of pulse waveforms must be accomplished 
about this reference point. 
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Figure 5: Sketch of EKG Pulse Form 

However, when the R pulse is detected, theP and Q 
pulses have already passed by. Thus, adequate 
memory of the P and Q pulses must be available so 
that they will be represented in the averaging which 
must be performed relative to the reference point 
on the R pulse, if detected. 

Since the HYDAC System operates with serial mem- 
ory, the speed necessary for the averaging must be 
obtained by adding the incoming period to the ac- 
cumulated sum of all previous periods in one rota- 
tion of the serial memory line. This demands that 
corresponding numbers to be added must be present 
simultaneously at the output of the sum memory and 
the input buffer memory. This, in turn, requires an 
alignment of the data in the two memories which, 
at worst, can take 256 memory cycles or about 500 
milliseconds. Thus, the arithmetic calculations take 
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Figure 6: Cardiogram Timing Diagrams 



500 ms or less to perform. Since cardiac periods 
lasting nominally from 500 ms to 2 seconds must be 
accommodated and, since it is not clear a priori how 
long a given record's period will be, a total of three 
input buffers are required. (See Figure 6 for timing 
diagrams.) These input buffer memories must be 
adequately controlled by a set of flip-flop state 
controls so that the appropriate memory either is 
being filled with data or is dumping data into a sum- 
mer leading to an accumulation or sum memory. 
Further, appropriate down-counting must be avail- 
able to detect an appropriate time interval following 
the R piilse. When this time, Tmin (the time of the 
shortest expected remaining period after the R 
pulse), is reached, a new memory begins accepting 



data. Loading of the original memory continues 
to time, Tjj^ax (the maximum time after the R 
pulse when significant parts of the waveform may 
still occur). This process is continued for n times 
(where n =16, 32 or 64) when the sum memory con- 
tains the desired average waveform. In addition to 
the EKG signal, the mean heart rate, T-rR > mean 
square heart rate, Tp^^^^, and mean square QRS 
time must be computed since these quantities are 
not available on the average waveform. One-word 
biiffers store the three associated pieces of input 
data to reduce access time to the 256 word buffers. 
The results then are punched out onto the paper 
tape as a more representative average waveform. 
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CONCLUSIONS 

A general purpose HYDAC program for the analysis 
of electrocardiogram data has been described. The 
major features of this program include: 

# The extension of this high-speed capability 
to permit on-line use of the computer pro- 
gram in complete system simulations that 
must proceed in real-time. 

# The capability for processing large quantities 
of data at high speed to provide statistically- 
significant results for meaningful cardio- 
graph analyses. 



9 The versatility provided by this general pur- 
pose approach... the output of Program I can 
be used as input to Program II or as input to 
any of a variety of other analytical devices, or 
the entire program can be modified for use in 
the analysis of other physiological processes. 



# The availability of adequate memory and 
timing controls to accomodate varying car- 
diac periods lasting nominally from 500 ms 
to 2 seconds and thus to account for individual 
differences. 



APPENDIX I 
THE MEASURED DATA 



Times... All times are measured from the beginning 
of the pulse whose letter is given first until the end 
of the pulse whose letter is given last. Thus, Tq^q 
lasts from the beginning of the Q pulse to the end 
of the S pulse. One exception to this definition is 
the timing signal Tj^r which equals the time between 
succeeding R pulses. 



for hST practically identical to the normal arith- 
metic average hei^t since the ST segment is gen- 
erally less than about 0.3 seconds. 

Variances... The variances of the Tq^gtime and the 
Tri^ time over a record are calculated using the 
relationship 



Amplitudes... The peak values of the P, Q, R, S and 

T pulses are measured if detected. Otherwise 
appropriate abandon search records are generated 
and recorded onto the paper tape. 

Height of the ST Segment. .The average height of the 
segment of the cardiogram between the S and the T 
pulses is computed on an integrator. The exponen- 
tially weighted average is used with an appropriate 
RC time constant ( 1 sec.). This generates a value 



2 2 

T = T - (T ) 
RR RR RR 



where = variance 



RR 



RR 



= mean or average of T 



RR 



= mean or average of T 



RR 
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APPENDIX II 



CRUCIAL EVENTS DETECTED ON THE AVERAGED WAVEFORM 
BY THE DECISION LOGIC 



Baseline Definition. ..The height of the EKGSD^y 

averaged waveform when its slope first lies 
between ±€ (a small number). 

Pulse Detection Definitions... (All events are expected 
in the following sequence): 

P Pu/se... Start =EKG3D > Vq (a small positive 
threshold). This is only recogniz- 
ed after the zero level baseline 
section has been detected and until 
'abandon point #1' is passed. 

Peak = ^ (EKG3D)^y =0 



d(EKG3D) 
R Pulse... Peak = = 0 



End =EKG3D<P, 



R negative slope = critical level - defines 
basic reference point. 



S Po/se... Start = EKG3Dav-Nc (Following the R 
pulse) 



End =EKG3Dav^Pc 



Peak =l(™)AV = o 



End = EKG3D>N, 



Q Pu/se.. .The search for the Q pulse begins after 
passing 'abandon point #1' or directly after 
P pulse ends, if it is detected. 



Start 



Peak = 



EKG3D^y<N(2 (a small negative 
threshold), and 



dt (EKG3D)^y<0 

d(EKG3D)AV ^ 0 
dt 



End = EKG3D>Ng 

Abandon search (A2) when EKG3Dav>Pc' 

Also start QT and QRS time measurements 
if they have not already begun by a de- 
tected Q pulse. 



Abandon search if S pulse does not occur 
within 160 ms past the R pulse. 



ST Segment... Lasts from the end of the S piilse or 
160 ms past the R pulse if the S pulse 
does not occur, until the beginning of the 
T pulse. 



T Pu/se... Start = EKG3Dav^Pc (After R pulse) 



dt 

End = EKG3D <P„ 



EAI 



ELECTRONIC ASSOCIATES, INC. West Long Branch, New Jersey 



ADVANCED SYSTEMS ANALYSIS AND COMPUTATION SERVICES/ANALOG COMPUTERS/HYBRID ANALOG-DIGITAL COMPUTATION EQUIPMENT/SIMULATION SYSTEMS/ 
SCIENTIFIC AND LABORATORY INSTRUMENTS/INDUSTRIAL PROCESS CONTROL SYSTEMS/PHOTOGRAMMETRIC EQUIPMENT/RANGE INSTRUMENTATION SYSTEMSAEST 
AND CHECK-OUT SYSTEMS/MIUTARY AND INDUSTRIAL RESEARCH AND DEVELOPMENT SERVICES/FIELD ENGINEERING AND EQUIPMENT AAAINTENANCE SERVICES. 



